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p I- ABSTRACT 

■ Context. Measuring the surface abundances of AGB stars is an important tool for studying the effects of nucleosynthesis and mixing 
' in the interior of low- to intermediate mass stars during their final evolutionary phases. The atmospheres of AGB stars can be strongly 
, aff'ected by stellar pulsation and the development of a stellar wind, though, and the abundance determination of these objects should 

therefore be based on dynamic model atmospheres. 

r~n ■ Aims. We investigate the eff'ects of stellar pulsation and mass loss on the appearance of selected spectral features (line profiles, line 

' intensities) and on the derived elemental abundances by performing a systematic comparison of hydrostatic and dynamic model 

, atmospheres. 

• . Methods. High-resolution synthetic spectra in the near infrared range were calculated based on two dynamic model atmospheres (at 

■ various phases during the pulsation cycle) as well as a grid of hydrostatic COMARCS models with eff'ective temperatures Teff and 
Mh' surface gravities log g over an adequate range. Equivalent widths of a selection of atomic and molecular lines (Fe, OH, CO) were 
Q \ derived in both cases and compared with each other. 

^ . Results. In the case of the dynamic models, the equivalent widths of all investigated features vary over the pulsation cycle. A consistent 

^ ■ reproduction of the derived variations with a set of hydrostatic models is not possible, but several individual phases and spectral 

^ ' features can be reproduced well with the help of specific hydrostatic atmospheric models. In addition, we show that the variations 

I in equivalent width that we found on the basis of the adopted state-of-the-art dynamic model atmospheres agree qualitatively with 
observational results for the Mira R Cas over its light cycle. 

H ■ Conclusions. The findings of our modelling form a starting point to deal with the problem of abundance determination in strongly 

^ ' dynamic AGB stars (i.e., long-period variables). Our results illustrate that some quantities such as the C/0 ratio can probably still 

^-H \ be determined to a reasonable accuracy, but the measurement of other quantities will be hampered by the dynamics. The qualitative 

OO . agreement with observations of R Cas opens promising possibilities for a forthcoming quantitative comparison of our synthetic spectra 

■ with observed ones of AGB variables in the globular cluster 47 Tuc. 

CO ■ 

• ■ Key words. Stars: late-type - Stars: AGB and post-AGB - Stars: atmospheres - Stars: abundances - Line: profiles Stars: variables: 
' general 

o : 
o. 

• • , 1, Introduction occurs concerning the spectral analysis. Stars in this part of the 

^ ■ HRD exhibit large amplitude radial pulsations. Being members 

Measurements of abundances in atmospheres of highly evolved ^^e variability class of so-called long period variables (LPVs), 

^ red giants contribute to the understanding of two major astro- jj^^^g objects change their visual brightnesses by up to several 

^ . physical aspects, namely the processes inside red giant stars such magnitudes over a few hundred days. The pulsating stellar in- 

. - - as nucleosynthesis and mixing, and, secondly, the enrichment of jg^j^rs lead to complex velocity fields in the stellar atmospheres 

the interstellar medium by freshly produced elements due to the Nowotny et al. 2010) and subsequently to phase-dependent 

mass loss of these stars. The determination of accurate element ^^^^^^ variations (asymmetric and split lines or emission 

abundances is complicated by the large number of spectral lines features) as seen in observed spectra of these stars (e.g. Hinkle 

from neutral atoms and molecules that give rise to heavy blend- ^^^2). Thus, a hydrostatic description of these objects is in- 

ing, especially in the optical range of the spectrum. The problem appropriate. However, it is this evolutionary phase where major 

IS less severe in the near infrared range, which has become in- nucleosynthesis and mixing processes take place, and therefore 

creasingly important for this kind of analysis. g^e^j jj^je^e^t in measuring stellar abundances exists. 

When a low or intermediate mass star evolves to the 

Asymptotic Giant Branch (AGB) phase, a second major problem In the past, a few attempts have been made to investigate 

the complicated profiles of individual lines observed in spectra 

* Operated by the Association of Universities for Research in of variable AGB stars on the basis of model calculations. Using 

Astronomy, Inc. under cooperative agreement with the National Science the dynamic models for pulsating red giants without mass loss 

Foundation of Bessell &. Scholz (il989) . Scholz (1992J showed how veloc- 
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ity fields in Mira photospheres distort line profiles of diff'erent 
species and aff'ect measurements of equivalent widths and curves 
of growth for abundance analyses from these spectra. This fun- 
damental paper, also referencing earlier works in this field, in- 
vestigated the possibility to "predict accurately the equivalent 
width of lines which are little aff'ected by outflow/infall veloci- 
ties by means of a conventional hydrostatic model in place of a 
dynamical model" (Scholz 1992). He found that lines formed in 
deep layers below the arising shock fronts may be accessible to 
a classical curve of growth study, but that in some cases this ap- 
proach seriously fails, so that abundance determination for large 
amplitude red variables remains a difficult problem. 

In this context, a major challenge for dynamic model atmo- 
spheres was to reproduce the line profile variations observable 
in spectra of Mira type variables. First results concerning this 
were presented by Bessell et al. (1996), who calculated syn- 
thetic line profiles (CO, Fe I) based on the models of Bessell 
& Scholz (1989). They could qualitatively reproduce some of 
the observed line profile characteristics (asymmetries, shifted, 
or doubled lines). 

McSaveney et al. (l2QQ7b attempted to derive abundances 
from high resolution spectra of pulsating AGE stars in the 
Magellanic Clouds by using the same type of model atmo- 
spheres. They succeeded in fitting the spectra of some stars 
at selected phases during the pulsation cycle. However, they 
had to restrict themselves to cases of relatively smooth atmo- 
spheric velocity fields. They were unable to derive abundances 
for model phases with complex atmospheric structures (shock 
waves), which cause distortions in the line profiles (emission 
components, line splitting). 

While the Australia-Heidelber g mod el atmospheres (e.g., 
Bessell & Scholz 1989. Bessell et al. [T996l Hofmann et aL iTWll) 
represent pulsating red giants without mas s loss, the Be rlin mod- 
els (e.g. Fleischer et al.T992, Winters et al. [T997t t|2000^ describe 
the optically thick dusty outflows of evolved C-rich AGB stars 
with h igh mass-loss rates. Based on the latter. Winters et al. 
(l2QQQb presented synthetic fundamental and first overtone CO 
line profiles. By comparing reasonably with observed ones for 
the obscured C-rich Mira IRC-i-10216, these spectral lines can 
be used to investigate the radial structure and the temporal vari- 
ation in the dusty envelope of these objects. 

The non-grey dynamic model atmospheres of Hofner et al. 
(I2QQ3D represent a combination of both approaches and describe 
the outer layers of an AGB star from the pulsating deep pho- 
tosphere (piston driven), to the dust forming region where the 
stellar wind is triggered and beyond to the region of the steady 
outflow (Nowotny et al. 12005 at . Keeping in mind the work of 
Scholz (1992), we conducted a similar exercise to investigate 
the problem of deriving abundances for LPVs by utilising the 
advanced dynamic models of Hofner et al. (12003b . 

The major aim of the work presented here was to study the 
efl'ect of atmospheric dynamics (pulsation, stellar wind) on the 
intensity of spectral lines for typical AGB stars. For this purpose, 
we compared selected spectral features in synthetic spectra com- 
puted on the basis of the aforementioned dynamic atmospheric 
models (Sect.|4l) with the corresponding results based on hydro- 
static model atmospheres (Sect.O. In a second step, we investi- 
gate observed high-resolution spectra of a typical Mira to see if 
the efl'ects we found by the modelling are also recognisable in 
the observational data (Sect.O. 



Table 1. Characteristics of the dynamic atmospheric models 
used for the modelling (P = pulsating, PM = pulsating and mass- 
losing). The notation was adopted from previous papers (Hofner 
et al. 120031 Nowotny et al. [2QTQl) . 



Model: 


P 


PM 


luminosity [L©] 


4000 


4000 


mass [M©] 


0.6 


0.6 


effective temperature [K] 


3500 


3500 


metallicity [Fe/H] [dex] 


-0.7 


-0.7 


C/0 ratio 


0.48 


0.48 


radius [Rq\ 


172 


172 


surface gravity log(g^ [cm/s^]) 


-0.25 


-0.25 


piston period P [d] 


200 


200 


piston velocity amplitude Awp [km s~^] 


2 


3 


luminosity amplitude parameter /l [1] 


4 


4 


bolometric variation Amboi [mag] 


0.86 


1.36 


dust mass absorption coefficient k^^^ [cm^ g~^] 




5 


mean mass loss rate (M> [Moyr~^] 
mean outflow velocity {u) [kms~^] 


1.9-10-^ 
19.5 



Table 2. Effective temperatures and surface gravities (ranges, 
stepwidths) of the grid of hydrostatic atmospheric models 
(COMARCS) used. 



C/0 


= 0.48 


TefF [K] 


2800. 


.4000 


A=50 






log(g [cm/s^]) 


-0.50 . 


. +0.50 


A=0.25 


C/0 


= 0.25 


Teff [K] 


2800. 


.4000 


A=50 






log(g [cm/s^]) 


-0.25 . 


. +0.75 


A=0.25 



2. Methods 

The starting point for the calculations are atmospheric structures, 
either from hydrostatic or dynamic models. Based on this, we 
calculate synthetic spectra and derive equivalent widths for a se- 
lection of atomic and molecular features. In our study, we fol- 
low the approach of an observational abundance study: to have 
a starting point for fitting an observed spectrum with a synthetic 
one, basic quantities such as effective temperature are derived 
from observed colours and subsequently the observed target is 
compared with the corresponding model. In our case, we also 
use a typical colour, (J-K), as an indicator of Tqq. That is, we in- 
vestigate how line intensities deviate in the dynamical case from 
the hydrostatic model, both related by a given colour. In doing 
so, we assume that (J-K) correlates with Tqq, keeping in mind 
that the term effective temperature may become questionable for 
the very evolved AGB stars featuring stro ng puls ations and mass 
loss (Baschek et al. lTWTl Nowotny et al. l20Q5bb . 

2.1. Model atmospheres 

The guideline for the parameters of our model atmospheres was 
a typical globular cluster LPV, namely the variable V3 in the well 
studied cluster 47 Tuc. This metal poor, short period mira was 
observed in the past by different authors (cf. Lebzelter & Wood 
[20051 Lebzelter et al. [2QQ5l t l2QQ6l and references therein)El The 
reason for choosing this object (and the corresponding models) 
was (i) the availability of spectroscopic data0 and (ii) the param- 



^ Approximate values found for V3 are: 
L = 4592 Lo, Teff = 3565 K, P = 192 d, M = 10-^-10-^ Moyr"^ . 
^ To be presented in a follow-up paper. 
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eters characterizing the star (comparably high Tq^, lower metal- 
licity) to reduce the confusing effect of line blending. 

For our simulations, we used the dynamic model atmo- 
spheres of Hofner et al. (2003 ). These represent a combined and 
self-consistent solution of hydrodynamics, frequency-dependent 
radiative transfei0 and - at least in the carbon-rich case - a de- 
tailed time-dependent treatment of dust formation and evolution. 
Being well suited to describing pulsating, mass-losing red gi- 
ants, these models were able to reproduce observational results 
such as low-resolution spectra of C-rich stars (Gautschy-Loidl 
et al. 120041) , line profile variations in high- resolution IR spec- 
tra (Nowotny et al. 2005a, 2005b||2010|), or interferometric data 
(Paladini et al. 2009 ) reasonably well. The parameters of the at- 
mospheric models (see Table [T]) were chosen to resemble our 
reference globular cluster mira. Starting with the same hydro- 
static initial model, we ended up with two diff'erent dynamic 
models - P and PM - due to the altered piston velocity amplitude 
Aup of the varying inner boundary (cf. Eqs. 1 and 2 in Nowotny 
et al. 120101) . While the first one represents a pulsating dust- free 
atmosphere, we encounter in the second case the development 
of a stellar wind on top of the pulsating photosphere. Plots of 
moving mass shells and the corresponding atmospheric struc- 
tures during a pulsation cycle for similar models of both types 
can be found in Hofner et al. (120031) or Nowotny et al. (12010 ). 
We note that in the case of oxygen-rich chemistry the physics of 
dust formation and m ass lo ss of AGB stars is not yet fully un- 
derstood (e.g. WoitkeEQOSl Hofner |2002]& [2008]). For M-type 
models such as those used here0 the dust-driven wind is there- 
fore induced by a simple parametrised description of the dust 
opacity first introduced by Bowen ("1988') and described in de- 
tail in Hofner et al. (2003; Eq. 19). We applied a maximum dust 
mass absorption coefl&cient of a:^^^ = 5 cm^ g"^ to generate such 
an artificial wind for model PM|3 

Since we wish to study the influence of atmospheric dynam- 
ics on spectral lines, we also computed a reference grid of hy- 
drostatic model atmospheres. We thereby followed the approach 
described in detail by Aringer et al. (120091) . Being guided again 
by 47 TUC-V3, a grid of COMARCS models was calculated 
with the aim of covering the full colour range of the dynamic 
models, as illustrated in Fig s. [2] and O The parameter ranges of 
the two subgrids with difl'erent C/0 ratios are given in Tabled 
Other model parameters were kept constant according to values 
found in the lit erature for 47Tuc-V3: =0.6Mq (Lebzelter 
& Wood 120051) . [Fe/H] = -0.7 d ex (C arretta & Gratton .1997 ), 
^ = 2.5kms"^ (Nowotny et al. 12010b . For the isotopic ratio 
^^C/^^C, we chose a typical AGB value of 30 (cf. Lebzelter et 
al.[2Q08l). 



^ The radiative transfer is solved at 64 frequency points and based on 
opacity sampling data of molecular opacities. 

^ For the C/0 ratios, the solar value of 0.48 was assumed, dynamic 
models for other C/0 ratios < 1 are currently under development. 

^ The critical value (i.e. corresponding to F = <3^adV<3grav = 1) of /Cd for 
the chosen set of stellar parameters (Table [B is 2 cm^ and the cho- 
sen K^^^ may appear to be a little high. Still, this is probably not an 
unrealistic assumption and comparable to what can be estimated for sil- 
icates in the context of dust-driven winds (see, e.g., Hofner 2008 ). Test 
calculations of another dynamic model with a reduced k^^^ of 3 cm^ g~^ 
(but higher piston velocity amplitude) led to a quite similar behaviour of 
an outflow strongly driven by radiative pressure. Moreover, not only is 
the dust absorption relevant for the wind, but also the effect of pulsation 
- simulated by the piston - providing kinetic energy input and lead- 
ing to enhanced densities in the cool outer layers. It should indeed be 
possible to obtain a stellar wind with an even smaller k^^^ for a higher 



2.2. Synthetic spectra and photometry 

For the spectral synthesis, we followed an approach that had 
already been applied very succes sfully i n the past, for details 
we refer to e.g. Nowotny et a l. (12005 af tJ20 10"), Aringer et al. 
(2002+2009), Lebzelter et al. (l2QQ8l) . or Gautschy-Loidl et al. 
(l2004D. 

Based on a given radial atmospheric structure (temperature- 
pressure) and the assumptions of chemical equilibrium as well as 
conditions of LTE, opacities for all necessary sources were cal- 
culated by using the COMA code (described in detail by Aringer 
2000, Aringer et al. |2009l Lederer & Aringer |2QQ9l) . Element 
abundances for the spectral synthesis were used consistent with 
the hydrodynamic models. We adopted the values for solar com- 
position from Anders & Grevesse (119891) , except for C, N, and 
O where we took the data from Grevesse & Sauval ("1994"). To 
obtain the correct metallicity (cf. Table[T]), the abundances of all 
elements heavier than H and He were then scaled by the same 
factor of 0.2 (as [Fe/H] = -0.7 dex corresponds to Z/Z© = 0.2). 
The C/0 ratios were subsequently adapted by increasing the 
amount of carbon (without changing the rest of the composi- 
tion). Following the descriptions given in Nowotny et al. (120101) , 
Doppler profiles were assumed for molecular lines and the cor- 
responding absorption coefficients were computed with the help 
of line lists by Goorvitch & Chackerian (119941) for CO and 
by Rothman et al. (2005; HITRAN database) for OH. For the 
atomic Fe lines, we followed the method of Lederer & Aringer 
(2009 ) and adopted full Voigt profiles with line data taken from 
the VALD database (Kupka et al.'2000'). All sources of the con- 
tinuum opacity listed in Lederer & Aringer ( 2009 ) were taken 
into account. In addition, we included absorption of H2O where 
necessary. Spectral features of water are present across the whole 
near-IR range and can severely influence the spectra for low 
effective temperatures and mass-losing objects as illustrated in 
Fig. [T] However, most of the existing H2O line lists show con- 
siderable differences in the positions of lines (Aringer et al. 
2008 ) when compared with observed spectra. Therefore, we ac- 
counted for water in a pseudo-continuous way in order not to in- 
troduce uncertainties related to the studied line profiles (blend- 
ing). The water absorption was calculated by using the (low- 
resolution) OS data from the SCAN data base (J0rgensen | 1 9971) , 
a detailed description of which can be found in J0rgensen et al. 
(|1993f t i2QQT1) . Spline fits through minimum values of appropri- 
ately grouped data points were then applied to determine the 
contribution of H2O to the actual high-resolution wavelength 
grid. This method ensures that the opacities are not overesti- 
mated (cf. Fig. 2.5 in Nowotny 2005 ) and infers lower limits to 
the quasi-continuous contributions. 

On the basis of the opacities Ky(A,r) for all depth and wave- 
lengths points determined in this way, the radiative transfer is 
solved in spherical geometry. As we wish to investigate the in- 
fluence of the complex atmospheric velocity fields on spectral 
lines, a radiative transfer code was applied that takes into ac- 
count the effects of macroscopic gas velocities. This code fol- 
lows the algorithm described by Yorke (1988 ) and has been used 
successfully in the past (e.g. Nowotny et al. 2005a+2010). Our 
aim was to study the behaviour of the individual line profiles. 
Thus, we computed synthetic spectra containing the spectral fea- 
tures chosen for our abundance studies (TableO with a high res- 
olution of 7? = 300 000. Rebinning to 7? = 70 000 led to spectra 
with resolutions comparable to observations with state-of-the-art 
IR spectrographs (FTS, CRIRES, Phoenix). The latter were then 
used for measuring equivalent widths (Sect. l2.3l) . Sample spec- 



4 



T. Lebzelter et aL: Abundance analysis for long period variables 




Fig. 1. Overview of the spectral features 
studied. Synthetic spectra based on three 
different COMARCS models (hydrostatic, 
i.e., no velocity effects) with log g = 0.0 are 
shown, for Teff^SSOOK (black), 3000 K 
(green), and 2800 K (red), respectively. The 
grey shaded areas indicate the ranges used 
to measure equivalent widths (cf. Table[3l). 
The spectra were normalised to the model 
continuum. Note that for the measurements 
a local pseudo-continuum was chosen. Some 
other spectral features are also identified for 
orientation purposes. Note that for decreasing 
temperatures spectral features caused mainly 
by H2O become prominent. For the modelling 
in this work, we accounted for water with 
a pseudo-continuous approach (see text in 
Sect. l2.2l ) to avoid uncertainties caused by line 
lists of limited precision in wavelength (in 
contrast to the spectra shown here, where line 
list data of H2O were utilised). 



2.313 2.314 



2.315 



2.366 



2.367 



2.368 



Table 3. Spectral features investigated. See Fig.[T] for an illus- 
tration of the feature locations. Column 3 gives an identification 
for each feature, which is used throughout this paper (contain- 
ing the corresponding central wavelength). Columns 4+5 list the 
borders of the spectral ranges used for determining equivalent 
widths. 



Species 


transition 


ID 


'^start [A] 


^end [A] 


Fe 




Fe 15536 


15535.5 


15536.4 


Fe 




Fel5538 


15537.6 


15539.1 


OH 


3-1 Pie 5.5 


OH 15540 


15539.1 


15540.3 


OH 


blend 


OH 15570 


15568.8 


15571.0 


OH 


2-0 Pie 11.5 


OH 15573 


15572.2 


15574.7 


OH 


2-0 Pif 11.5 


OH 15576 


15575.7 


15577.0 


i^CO 


3-0 head 


CO 3-0 band head 


15581.0 


15585.7 


i^CO 


2-0 R19 


i^CO 23137 


23136.3 


23139.6 


i^CO 


2-0 R19 


i^CO 23360 


23658.8 


23662.5 



tra based on hydrostatic COMARCS models (thus no velocity 
effects) are shown in Fig.[T] 

In addition to these high-resolution spectra containing the 
line profiles, we also computed low-resolution spectra and syn- 
thetic photometry following the approach of Aringer et al. 
(I2QQ9I ). Opacity sampling spectra covering the near-IR with res- 
olutions of R = 10 000 were calculated on the basis of several at- 
mospheric models. By convolving these spectra with filter trans- 
mission curves and applying adequate zeropoints, we obtained 
synthetic photometry. The resulting {J-K) colours (subsequently 
used to correlate different models) for the dynamic model atmo- 
spheres during the pulsation cycle and for the grid of hydrostatic 
COMARCS models, are shown in Figs. [2]and[3l respectively. 



2.3. Measurement of equivalent widths 

Elemental abundances in the atmospheres of stars can be de- 
rived from the intensity of spectral features since the number 
of absorbers will directly affect the depth of a line (except for 
saturated lines) apart from effects of temperature and electron 
density. Since we wish to investigate the effect of atmospheric 
dynamics on the derived abundances, we have to look for a quan- 
tity describing the line intensity. The standard concept in the 
context of abundance measurement is the equivalent width of a 
line, which is related to the abundance of the parent element by 
the curve of growth (e.g. Castelli & Hack 1990). The spectra of 
cool giants with their severe blending of adjacent molecular lines 
would certainly favour the application of a full spectral compar- 
ison to study the abundances of the elements on the surface (see 
e.g. Lebzelter et al. 2008 ). However, as we see below, stellar pul- 
sation leads to considerable modifications of the line profile with 
the appearance of emission or doubled absorption components 
(cf. Nowotny et al. l2005al 12010b . Therefore, a complete fit of a 
spectrum of a mira-like variable in a particular spectral range 
with hydrostatic model spectra will fail anyway (at least during 
a good portion of the pulsation cycle). This makes a comparison 
of hydrostatic and dynamic models difficult. With the help of the 
equivalent width, we can at least measure the amount by which 
the total line intensity is changed by pulsation. 

We obtained equivalent widths for a selection of spectral 
features of various species. For this, we first imported the syn- 
thetic spectra into IRAK To allow for an automatic and objective 
measurement, we removed all velocity shifts from the dynami- 
cal spectra by shifting the investigated features to zero veloc- 
ity (while retaining the Doppler- modified line profiles). For this 
shift, we used the line core as a reference point, except in cases of 
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Fig. 2. (J-K) colour variation versus bolometric phase for the 
dynamic models with (PM, filled symbols) and without (P, open 
symbols) wind. 




4000 3800 3600 3400 3200 3000 2800 

Teff [K] 



For our study, we selected lines in the H- and ^-band. These 
wavelength ranges have played a key role in the study of pho- 
tospheres of AGB stars since the development of high resolu- 
tion near-infared spectrographs. Within this wavelength range, 
we chose a combination of strong and weak lines, and atomic 
and molecular features including a molecular band head for our 
analysis. Many of the selected lines had already been used in 
the abundance analysis of red giants (e.g. Smith et al. 120021 
Lebzelter et al. 12008b . For each feature we measured the same 
area around the central wavelength in all spectra. The borders 
are listed in TableO Examples illustrating the measured range 
are shown in Fig.[T] 

A few uncertainties in this kind of measurement have to be 
considered. First, variable line asymmetries as well as a lim- 
ited number of points in the pseudo-continuum can introduce 
small inconsistencies in the handling from one spectrum to an- 
other. We ran some tests and found possible deviations below 
2% of the measured equivalent width. A second and more severe 
problem is blending, although we tried to avoid strongly blended 
lines in our selection of spectral features. The problem emerges 
mainly for cool stars (Tqq <3000 K) with the appearance of H2O 
and CN lines in the spectrum. We return to this point below. For 
the dynamic models, the occurrence of line doubling and emis- 
sion components is a problem. The former can lead to signifi- 
cant broadening of the feature requiring us to measure a wider 
range about the central wavelength. This, however, leads to an 
inclusion of any (typically blueshifted) emission component re- 
ducing the equivalent width of the feature at other phases. When 
measuring a line in an observed stellar spectrum, one may delib- 
erately exclude this emission component provided the spectral 
resolution is high enough. This would lead to qualitative and 
phase-dependent diff'erences between our model results and ob- 
servations. Emission components close to the line core do not re- 
quire special consideration since any measurement in observed 
spectra would not be able to distinguish them either. 



Fig. 3. Relation between T^ff of the hydrostatic COMARCS 
model atmospheres and the resulting (J-K) colours. Diff'erent 
line styles show a variation in logg from light grey dotted 
(log^=-h0.5), light grey dashed (-hO.25), solid black (0.0), dark 
grey dashed (-0.25), to dark grey dotted (-0.50). In the online 
colour version, this corresponds to cyan (-hO.5), green (-hO.25), 
black (0.0), red (-0.25), and blue (-0.5), respectively. The in- 
sert shows the same relation for log g =0.0 but two diff'erent C/0 
ratios, illustrating the influence of a change in this parameter. 



line doubling, where we used instead a mean value between the 
two components. For the shifted spectra, we afterwards verified 
that the complete feature falls within our range of measurement 
and corrected the shift if necessary. All spectra were normalized 
to the local pseudo-continuum instead of the model continuuni^ 
to resemble an observer's approach as realistically as possible. 
The pseudo-continuum was defined by placing a straight line 
through the highest points in the spectrum excluding obvious 
emission peaks (as present in the case of the dynamic models). 
The equivalent widths were then measured using the IRAF tasks 
splot and bplot. 



^ In contrast to the pseudo-continuum approach usually applied by 
observational studies, modelling offers the possibility to normalize the 
synthetic spectra relative to a computation for which only the continu- 
ous opacity was taken into account {F IF com', cf. Fig.[B- 



3. Dependence of hydrostatic line sensitivities on 
stellar parameters 

We describe the dependency of the intensity of various spectral 
features on (i) eff'ective temperature (reff), (ii) surface gravity 
(logg), (iii) C/0 ratio, and (iv) ^^C/^^C ratio as derived from 
synthetic spectra based on hydrostatic models. The results are 
grouped by the element or molecule producing the feature. 

As can be seen from Fig.O the roughly linear relation be- 
tween Teff and {J-K) stops at a temperature of 3200 K. For 
lower temperatures, the near-infrared colour becomes almost 
constant and even shows a slight linear trend towards lower 
values. This phenomenon has also been found for carbon rich 
stars by Aringer et al. ( 2009 ) and is caused by the behaviour of 
diff'erent molecular features in the region of the filters. For the 
same reason - but of course due to diff'erent molecular species 
(mainly H2O) - a similar trend is seen in the oxygen-rich models. 
However, AGB stars in this low-temperature range are subject to 
mass-loss processes and the eff'ect induced on the colour by the 
subsequent circumstellar reddening may be much larger than the 
one caused by the aforementioned behaviour of molecular fea- 
tures. This is especially true for C-rich stars and high mass-loss 
rates (Aringer et al. 12009b . 

As shown below, the intensities of the spectral features also 
continue to increase below 3200 K (except for the iron lines, 
which decrease in strength with decreasing Teff). This results in 
a turnaround in the equivalent width curves at the high {J-K) 
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Fig. 4. Equivalent widths in [A] of the Fe line at 15 536 A mea- 
sured from synthetic spectra based on hydrostatic model atmo- 
spheres against the corresponding (J-K) values. Same line style 
coding for various logg values as in Fig.[3l Details of the range 
in wavelength used around the line centre are listed in Table[3] 
and illustrated in Fig.[T] 

values. We decided to plot the relations versus the near infared 
colour for two reasons. First, the (J-K) colour is a quantity that 
is directly measurable. It is thus easier to compare the modelling 
results with observations (see below). Second, the comparison 
with dynamic models (performed in the next section) requires 
the use of a global parameter such as the colour, because the 
concept of an effective temperature is questionable for extended 
dynamical structures. 

3.1. Fe lines 

Two iron lines were selected for the analysis to investigate the 
behaviour of lines typically used to derive stellar metallicities 
(e.g. Smith et al. 200 2). The two lines originate in the same ex- 
citation level (5.642 eV). Accordingly, both lines show a similar 
trend with Tqq and log g, there being a clear dependence of the 
equivalent widths on both parameters as illustrated for one of 
the lines in Fig.lH The line at 15 536 A is only mildly affected 
by blending, although a CN feature becomes strong for the low- 
est temperatures studied. The second Fe line close to 15 538 A, 
which behaves qualitatively identically, is blended with a nearby 
CN line producing the asymmetric line profile clearly visble in 
Fig.[T] A neighbouring OH line mildly affects the red wing of 
the line profile. As can be seen from Fig.[T] a large fraction of 
the observed decrease in the line strengths at low temperatures 
is caused by a depression of the overall continuum. 

3.2. OH lines 

Lines of the OH molecule produce very prominent features in 
the //-band spectra of cool giants. We selected two lines from 
the 2-0 transition and one line from the 3-1 transition for our 
study. In addition, we included a blend of OH lines across the 
same wavelength region in our analysis. A CN line is located 
in the blue wing of this blend. The blue end for measuring the 
equivalent width was defined to include the complete blend also 
in the dynamical case, but this means that most of the CN line is 
also covered (see Fig.[T]). 

All these lines behave identically. They show a strong depen- 
dency on (J-K), i.e. Teff , while the change in the line intensity 
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Fig. 5. Equivalent widths in [A] of the OH blend at 15 570 A 
measured from synthetic spectra based on hydrostatic model at- 
mospheres against the corresponding (J-K) values. Same line 
style coding for various log g values as in Fig. [3] The dash-dotted 
line indicates the relation for C/O=0.25 (standard value: 0.48). 




0.9 1 1.1 1.2 1.3 

J-K [mag] 

Fig. 6. Same as Fig.[5]for the OH line at 15 576 A. 



with log g is comparably small. This is shown for two of the lines 
in Figs. [Hand [6] The OH blend (Fig.O including also a CN line 
does not behave differently, the feature being dominated by OH. 
We also derived equivalent widths for a model with log g = 0.0 
and C/0 = 0.25 to study the effect of a change in the C/0 ratio. 
The other parameters were left unchanged. The result is plotted 
as a dash-dotted line in the figures. The effect is visible for the 
OH 2-0 line at 15 576 A, but it is rather small. In this context, 
we note that for our modelling - simulating the change caused 
by third dredge up - the C/0 ratio is changed by increasing or 
decreasing the carbon abundance, while the oxygen abundance 
remains constant. 

3.3. CO features 

Features of CO with its well expressed band structure can be 
easily identified in the H- and /T-band spectra of stars of spec- 
tral type K or later. They are frequently used to derive stellar 
parameters and abundances or to study the atmospheric dynam- 
ics (e.g. Hinkle et al. 1982, Lebzelter [T9991). The three CO fea- 
tures selected here represent different spectral signatures of this 
molecule: the (3-0) band head of ^^CO located in the //-band; 
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Fig. 7. Equivalent widths in [A] of the ^^CO 3-0 band head (mea- 
sured between 15 581 A and 15 585.7 A) measured from syn- 
thetic spectra based on hydrostatic model atmospheres against 
the corresponding {J-K) values. See Fig.[T] for the wavelength 
range covered. Same line style coding for various \ogg val- 
ues as in Fig.O The dash-dotted line indicates the relation for 
C/O=0.25 (standard value: 0.48). 
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Fig. 8. Equivalent widths in [A] of the ^^CO line at 23 137 A 
measured from synthetic spectra based on hydrostatic model at- 
mospheres versus the corresponding {J-K) values. Same line 
style coding for various log g values as in Fig. [3] The dash-dotted 
line indicates the relation for ^^C/^^C=10 (standard value: 30). 
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Fig. 9. Same as Fig.[8]for the ^^CO line at 23 660 A. 



a strong single ^^CO line in the ^-band; and finally a line from 
^^CO in the same wavelength range. 

As in all other cases, the equivalent widths of the three fea- 
tures were derived for various Teff and logg (Figs.|7] to [9]). In 
addition, the CO band head was also measured for a model with 
logg =0.0 and reduced C/O=0.25. To investigate the eff'ect of a 
change in the carbon isotopic ratio, we also calculated synthetic 
spectra with ^^c/ 13C=10 instead of 30 (log g =0.0). 

The strength of the CO band head (Fig. [7]) depends on all 
three parameters Teff , log and C/0. As one would expect from 
simple considerations concerning molecule formation, for the 
star investigated here the intensity of the C/0 bands increases 
with the C/0 ratio. Quite a diff'erent behaviour can be seen for 
the ^^CO line at 23137 A. The equivalent width is almost con- 
stant for a wide range of {J-K), but below 3400 K a strong de- 
pendency on Teff is observed. We note that, in contrast to the sit- 
uation around 1.5 yum, the line strengths of the features at 2.3 \im 
are not aff'ected by a change in the continuum level. The ^^CO 
line investigated is almost independent of logg. An exception 
is the model with the lowest surface gravity, which seems to be 
off'set from the other models. The reason for that behaviour is 
unclear, but a weakening of the lines due to sphericity eff'ects 
occurring in very extended atmospheres is a likely explanation 
(Aringer et al. 12009b . Finally, the ^^CO line again shows a trend 
of increasing line strenghts with T^^ and a mild dependency on 
the log g value (especially for log g <0.0, see Fig.O. The change 
with Teff seems to be intensified by a neighbouring line of CN, 
which appears at the lower end of the temperature range studied. 
A change in the carbon isotopic ratio has a strong impact on the 
intensity of the studied ^^CO line. 

4. Comparing dynamical and hydrostatic line 
intensities 

We summarize our findings for hydrostatic models as follows: 

- IMost features are sensitive to changes in T^^. The only ex- 
ception is the ^^CO line at 23 137 A, which becomes temper- 
ature sensitive for low temperatures only. The two studied 
lines of Fe are also less sensitive to Teff changes for higher 
temperatures. 

- The dependency on T^^ is a combination of a change in both 
the line intensity and the level of the pseudo-continuum due 
to H2O. 

- At the lowest temperatures, we observed CN line blending 
for most of the features studied. 

- A change in surface gravity has on average a smaller influ- 
ence on the feature strengths. The most pronounced changes 
in this parameter were found for the CO 3-0 band head and 
the iron lines. 

With this study of the hydrostatic case we defined a refer- 
ence grid with which to compare the dynamic models. We again 
discuss the effects of pulsation on a line by line basis. 

4A. Fe lines 

We first note here the overall behaviour of dynamic models in the 
{J-K) versus equivalent width plane, which is already clearly 
visible for the Fe lines in Figs.[TOl and [TTJ for a given {J-K) 
value, the line intensity is diff'erent for the falling and rising 
branch of the light curve. As a result, the model values loop 
through the diagram. This loop is more pronounced in the case 
of the model with mass loss {FM) than for the model without 
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Fig. 10. Change in the intensity of the Fe line at 15 536 A over 
the calculated pulsation cycles versus {J-K) for two dynamic 
models. The model with mass loss (PM) is shown with solid 
symbols, while the model without mass loss (P) is marked with 
open symbols. Data points are connected to illustrate the path 
of the change over the pulsation cycle. Two cycles are plotted 
for model PM and four for model P. The small arrow indicates 
the direction of the path and starts at the maximum bolometric 
phase (0^0/ =0). The solid line shows the hydrostatic model for 
log^ = -0.25 for comparison (see Sect.©. Compare with Fig.|4l 
for various values of log g. 



mass loss (P). Small changes are visible between the pulsation 
cycles calculated, but the general behaviour remains the same. 

In the dynamcial case, the {J-K) - equivalent width relations 
for the iron lines show a similar arc-like shape as in the hy- 
drostatic case. The maximum line intensity is reached at nei- 
ther maximum nor minimum light. However, the two iron lines 
clearly differ in their behaviour. Fe 15 538 A seems to vary in a 
way similar to that in the hydrostatic models. The variation in 
line intensity over the pulsation cycle for model PM may be de- 
scribed by a combination of a hydrostatic model of log g =0.0 
for the first half of the light curve and a hydrostatic model of 
somewhat lower surface gravity for the second half of the light 
variation. Model P is off'set to somewhat higher line strengths. 

The Fe 15536 A line (Fig.fTOl) shows a trend clearly diff'erent 
from the other iron line and from the hydrostatic reference case. 
Here the line strenghts are weaker than in the COMARCS model 
and would correspond to a higher surface gravity. In addition, the 
shape of the relation is diff'erent. Model P again exhibits stronger 
lines at a given (J-K) than model PM. For (J-K) >1.2, model P 
approaches the hydrostatic relation. 

4.2. OH lines 

All four OH lines investigated exhibit similar behaviour over the 
pulsation cycle. As for the iron lines, the model produces a loop 
in the diagram over one pulsation cycle. The dynamic model P 
again exhibits stronger lines than model PM. The diff'erence be- 
tween the two models changes from line to line. 

Both dynamic models have similar dependences on (J-K), 
i.e. Teff , as our hydrostatic reference models (Figs. [12] to [T4l) . 
During half of the pulsation cycle, we can see an off'set between 
model PM and the hydrostatic case. As the hydrostatic line in- 
tensity is almost independent of the surface gravity (cf. Figs. [5] 
and [6]), this off'set cannot be compensated for by choosing a dif- 
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Fig. 11. Same as Fig.[l0lfor the Fe line at 15 538 A. 
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Fig. 12. Same as Fig.[l0lfor the OH line at 15540 A. 



ferent logg value as in the case of the iron lines. To bring the 
dynamic model in alignment with a series of hydrostatic mod- 
els, one would be inveigled to change the C/0 ratio. However, 
the OH line at 15 540 A (Fig. [12]) shows a deviation from the hy- 
drostatic case towards the opposite side of the other OH lines, 
therefore requiring an inverse change of C/0. 

The diff'erent slope seems to be a result of a complex mix- 
ture of various components of the investigated lines, which are 
difficult or even impossible to distinguish in the measurement. 
One obvious feature is the appearance of an emission compo- 
nent between 0^^/ =0.26 and 0.41 (corresponding to the lower 
branch of the loop between (J-K) - 1.1 and 1.28; see Fig. [14]) 
on the blue wing of the OH 15573 and OH 15576 A lines. This 
emission component is not excluded from the measurement of 
the equivalent width and may be responsible for at least part of 
the depression of the line intensity measured. 

4.3. CO features 

The most remarkable loops in the diagram EW versus (J-K) 
corresponding to model PM are found for the features of CO. For 
the CO band head (Fig. [15]), model P (no outflow) and part of the 
pulsation cycle of model PM show a very similar relation of band 
strength with (J-K) as the hydrostatic models. Owing to the high 
sensitivity of the CO band head to the surface gravity (note the 
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Fig. 13. Same as Fig.[l0lfor the OH blend at 15570 A. Compare 
with Fig. 151 




J-K [mag] 



Fig. 14. Same as Fig.[TOlfor the OH line at 15576 A. Compare 
with Fig.O 



different scale on the y-axis in Fig. [7]), the whole loop made by 
model PM can be reproduced by a change in logg of 0.25. A 
modification of the C/0 ratio from 0.48 to 0.25 has an about 
three times larger eff'ect on the band strength than the changes 
during the pulsation cycle. 

For the two individual CO lines studied in the ^-band, the 
situation is completely diff'erent, especially for model PM (with 
outflow). For the ^^CO line at 23 137 A, we find only very short 
parts of the pulsation cycle where the derived line intensity is in 
a similar range as that found for the hydrostatic models. This de- 
viation from the hydrostatic case is caused by the occurrence of 
line doubling over most of the pulsation cycle. However, this 
does not lead to a simple broadening. As can be seen from 
Fig. [161 the line intensity at some phases is also weaker than in 
the hydrostatic case. This is caused by an emission component 
that is visible on the blue side of the line profile for almost half 
the light cycle. For illustration, we compare in Fig.[T71the corre- 
sponding line profile of this CO line derived from a hydrostatic 
model and from two phases of model PM at similar (J-K). The 
two chosen phases represent the two branches of the bolometric 
light curve. 

Model P resembles the hydrostatic case more closely. Line 
doubling is observed less often, but again a blue-shifted emis- 




J-K [mag] 



Fig. 15. Same as Fig.[TOlfor the ^^CO 3-0 band head. Compare 
with Fig. [71 

sion component can be seen. Here the diff'erences between hy- 
drostatic and dynamic models are greatest at minimum light. 

Line profile variations such as those havee frequently been 
reported in the litera ture (e.g. Hink le et al. 11982 1 Lebzelter et al. 
[T^ Nowotny et al. 12005 all2005bb . The low excitation 2-0 tran- 
sitions of ^^CO are strong lines that form over a large part of the 
stellar atmospheres and also in the outflow of AGB stars. This 
makes the difl'erence between line intensity variations in models 
P and PM understandable because the latter model has a larger 
extension and introduces an additional redshifted component to 
the line profile. 

For the ^^CO line studied, we also find a behaviour deviating 
from the hydrostatic case (Fig. [181). Model P shows a trend very 
similar to the hydrostatic models over a large fraction of the pul- 
sation cycle. The line intensity in model PM follows the hydro- 
static relation for about half the pulsation cycle, but for the other 
half the lines are much weaker. This cannot be accounted for by 
a change in log g alone and would require an enormous change 
in the isotopic ratio. The occurrence of an emission component 
in the line profile is responsible for this weakening as in the case 
of the ^^CO line. The emission is visible at several phases as a 
separate component, but is probably hidden by a general weak- 
ening of the line during a much larger part of the pulsation cycle. 
If we take out the emission component by artificially setting all 
values above the pseudo-continuum to the continuum level, the 
loop becomes less extended in the lower right corner of Fig.[T8l 
Contrary to the case of model P, a flattening of the loop cannot 
be achieved in this way. 

In Fig.[T8l we also see noticeable cycle-to-cycle difl'erences, 
especially in the case of model PM. This is an interesting re- 
sult keeping in mind that our model atmosphere is modulated 
by a strictly sinusoidal variation in the inner boundary. This be- 
haviour was previously noted for other output quantities (e.g. de- 
gree of dust condensation) from this kind of models (Nowotny 
et al. l2005a[ l2005b) . For models based on self-excitation pulsa- 
tion or real stars (see next section), even larger cycle-to-cycle 
difl'erences can be expected (Lebzelter et al. 120011) . 

5. Results from observations 

It is naturally it is a critical question whether the behaviour found 
for various lines in the dynamic models are also visible in obser- 
vations. For the hydrostatic case, a number of successful com- 
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Fig. 16. Same as Fig.[l0lfor the ^^CO line at 23137 A. Compare 
with FigM 
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Fig. 17. Comparison of the line profile of ^^CO 23 137 in a hy- 
drostatic model (Tqq =3650 K, log g =0.0) and two phases of the 
dynamic model PM. All three show a very similar (J-K) value 
close to 1.1. 
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Fig. 18. Same as Fig.[l0lfor the ^^CO line at 23 660 A. Compare 
with Fig.[9l 



parisons between observed and synthetic photometry and spec- 
tra exist (e.g. Loidl et al.l200U Aringer et al. 2002, Lebzelter et 
al. 12008b . For the dynamical case, it could also be shown that the 
models reproduce observations of C-rich stars reasonably well 
(e.g. Gautschy-Loidl et al. 2004, Nowotny et al. 2005b). 

Long time series of high resolution spectra of LPVs are rare. 
Here we use a time series of FTS spectra in the H- and ^-band 
of the mira RCas. This object is a typical mira variable in the 
solar neighbourhood with a period of 43 1 ^ and a visual bright- 
ness amplitude up to 9 magnitudes. At minimum light, the star 
reaches a spectral type of MIO (Lockwood 11970 , Wyckoff* & 
Wehinger J971i) . Its dusty envelope has been used to d erive a 
mass-loss rate of a few lO'^Moyr"^ (Schutte & Tielens [19891 
Weigelt & Yudin 2001). This agrees with an older mass-loss 
rate derived from thermal circumstellar CO emission (Loup et 
al. 11993b . The atmospheric kinematics of RCas was described 
by Hinkle et al. (1984"). These authors note that around mini- 
mum light the near infrared CO line profiles investigated by them 
seem to be strongly aff'ected by blending of the CO lines with 
H2O lines. Strong water absorption bands were also reported by 
Truong-Bach et al. ([T9991) and Aringer et al. (l2Q02b in the mid- 
infrared range. Clearly R Cas is described by parameters quite 
diff'erent from our models, and we wish to explicitly state that 
the dynamic model atmospheres (Table[T]) were not designed to 
resemble this specific object. We did not attempt to reproduce the 
spectrum of R Cas. The main goal of this study is to identify the 
basic problems and features of abundance analysis in dynamical 
atmospheres. Our comparison shall be made on a qualitative ba- 
sis only. The main reason for choosing this set of spectra is that 
to our knowledge it is the most extensive time series of high- 
resolution near-infrared spectra of an M-type mira available. 

Our time series was obtained by one of us (KH) between 
August 1983 and September 1986 with the FTS spectrometer at 
the Kitt Peak National Observatory 4 m Mayall telescope (Hall et 
al. n979i) . It covers two complete light cycles of this star. A stan- 
dard data analysis described extensively in previous papers (e.g. 
Hinkle et al. 1982, 1984 ) was applied to the spectra. Additional 
analysis was performed within IRAF using the same procedures 
as in the case of the synthetic spectra. 

Phases within the pulsation cycle were attributed to the in- 
dividual observations by using the AAVSO lightcurve (Fig. [2T1) . 
The figure clearly indicates that by chance our time series covers 
an untypical behaviour of the star with an outstandingly bright 
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Fig. 19. FTS spectrum of R Cas around the OH lines obtained at 
JD 2 445 575 (corresponding to 0v = 0.59). The line indicates the 
chosen location of the pseudo-continuum. Arrows mark the three 
OH lines at 15 570, 15 573, and 15 576 A, respectively (cf. upper 
right panel of Fig.[T]). Note the clear indication of line doubling 
in the two stronger OH lines. The CO 3-0 band head can be seen 
on the right side of the plot. 



maximum during the first cycle and an unusually weak max- 
imum in the second. The third maximum plotted represents the 
typical amplitude observed for this star. That this behaviour is re- 
flected in an untypical atmospheric dynamics was demonstrated 
by Lebzelter et al. (1200 l i This has to be kept in mind when 
discussing the line behaviour in R Cas. 

{J-K) values obtained parallel to our time series do not ex- 
ist. The atypical light change during that time does not allow us 
to use near-infrared light curves obtained during some other pul- 
sation cycle (e.g. Kerschbaum et al. 2001 , Nadzhip et al. 2001 ). 
The only quantity directly available is the visual brightness from 
the AAVSO light curve. Assuming that (/-^) changes approxi- 
mately in phase with the visual brightness, we use this indicator 
in our investigation. The V value was derived from a polynomial 
fit to all AAVSO data. No correction for the diff'erence between 
visual measurements and photometric V-band magnitudes was 
applied, but the eff'ect on the shape of the light curve should be 
very small. 

We attempted to measure spectral features in the R Cas data 
in a similar way (i.e. using the same width to determine the 
equivalent widths of lines) as in the synthetic spectra. We natu- 
rally had some difficulties in this case: the definition of the con- 
tinuum is far more complicated in the observations than in the 
synthetic spectra, where emission peaks can be easily identified 
from the nominal continuum level. We illustrate our definition of 
the continuum in Fig. [191 As for the synthetic spectra (Sect. l2.3l) , 
the observed spectra were shifted in wavelength (while preserv- 
ing the Doppler-modified profiles) to allow an automatic mea- 
surement of equivalent widths. For comparison, we chose fea- 
tures that are only slightly contaminated by telluric absorption 
and assume that the influence of the Earth's atmosphere is neg- 
ligible. Finally, our data are slightly afl'ected by the instrument 
profile for which we also did not apply any corrections. 

The most remarkable features in the case of the synthetic 
spectra are certainly the loops derived from the dynamic mod- 
els, especially those models that produce an outflow of material. 
For a given near-infrared colour, two line strengths are found 
depending on the star being on the rising or falling part of the 
light change. R Cas is obviously also a mira with a substantial 
mass loss. In Fig.[20l we show the behaviour of the OH line at 
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Fig. 20. Equivalent widths in [A] of the OH blend at 15570 A 
measured from our time series of FTS spectra of the mira R Cas 
(triangles) against the corresponding V values taken from the 
AAVSO light curve. The small arrow indicates the path from 
the first visual light maximum to minimum. For comparison, we 
also added the results from model P (open circles) and model PM 
(filled circles). Note that the model parameters are very diflferent 
from the parameters of R Cas. 

15 570 A versus visual brightness. Consecutive data points in the 
light cycle are connected by lines. Altering the continuum level 
by 2 % changes the measured equivalent width by about 0.03 A. 
Other spectral features checked by ourselves show a similar be- 
haviour. 

For an easier comparison, we also show the behaviour of the 
equivalent width with visual brightness derived from our mod- 
els P and PM. The absolute visual magnitudes predicted by the 
models were shifted to roughly coincide with the V-brightness of 
R Cas. It is clear that the observed behaviour in R Cas is in qual- 
itative agreement with that found in our study using synthetic 
spectra. However, we note that the loop over the light cycle goes 
in the inverse direction (see the arrow in Fig. [20|. At the mo- 
ment, the reason for this diff'erence is not clear. Figure [21] shows 
a strong correlation between the feature intensity and the visual 
brightness. The diff'erence between the first and the second ob- 
served light cycle mentioned above is therefore clearly reflected 
in a corresponding change in the equivalent width. As a conse- 
quence, the loops resulting from these two light cycles (Fig.l20l) 
have difl'erent extensions. 

The occurrence of loops in a digram of equivalent width ver- 
sus either {J-K) or V corresponds to a difl'erent relation between 
the two quantities on the rising and the falling branch of the light 
cycle. This can be seen in both Fig.[22lfor the model case and 
Fig.[2T]for the observations. 

6. Conclusions and outlook 

The aim of our study has been to explore the possibilities and 
caveats of deriving stellar abundances from observed spectra of 
AGB stars that are strongly afl'ected by pulsation and mass loss. 
We can summarise our findings as follows: 

- There is not a single phase, neither for model P nor for 
model PM, where all investigated features would be in agree- 
ment with a single hydrostatic model at the same time. This 
suggests that there is a fundamental problem in fitting the 
complete spectrum of a dynamical star with one hydrostatic 
model. 
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Fig. 21. Equivalent widths in [A] of the OH blend at 15570 A 
(large dots, left y-axis) and the AAVSO light curve (small dots, 
right y-axis) for R Cas against time expressed in JD. 




Fig. 22. Equivalent widths in [A] of the OH blend at 15570 A 
(solid symbols, left y-axis) and (J-K) colour (solid line, right 
y-axis) for model PM against the bolometric phase. Phase has 
been assigned to the light maximum 



For several features, the dynamic models, especially 
model PM, touch areas in the (J-K) versus equivalent width 
plane that seem to be unreachable by hydrostatic models 
even when modifying the chemical composition. Further 
tests are necessary to test the effects of e.g. metallicity on 
this finding. 

Spectral peculiarities due to stellar dynamics seem to be less 
significant in the //-band than in the /T-band. 
For dynamic atmospheres with mass loss (represented by 
model PM), the rising, post-minimum branch seems to be 
more promising for deriving element abundances with the 
help of hydrostatic models. This is, however, untrue for some 
of the features studied (namely the iron lines, OH 15540 and 
i^CO 23137). 

For dynamical atmospheres without mass loss (model P), the 
loops are less extended and the detected trends of equiva- 
lent width with colour seem to resemble the hydrostatic case 
in many respects. For several features, we find an off'set be- 
tween the relation for model P and the hydrostatic reference 
case with log g = -0.25. This off'set, however, is not the same 
for all features. 



- Including the emission lines in the measurement of the 
strength of a spectral feature leads to an extension of the 
loop completed in the colour versus equivalent width plane. 
For high resolution spectroscopy, it is therefore better to ex- 
clude any clearly separated emission component when de- 
riving the stellar abundance. If low resolution spectra were 
used, in which the emission component cannot be detected 
separately, we expect to see extensive loops such as those 
found in this study. 

- Our results are qualitatively similar to the findings of Scholz 
(1992) and others. 

We performed some tests to see if a combination of vari- 
ous spectral features (e.g. a ratio of the equivalent widths of two 
lines) could reduce the complexity of the variation with colour. 
But as the size and slope of the loops are diff'erent for each of 
the features observed, none of the combinations were found to 
be useful. 

Deriving stellar abundances for pulsating stars without an 
outflow (simulated in this study by model P with a bolomet- 
ric light amplitude of 0.86"^^^) may well be possible if we take 
into account a set of ofl'sets that still need to be determined. 
Abundance determinations for these kind of stars performed so 
far are probably afl'ected by systematic error, depending on the 
feature used. 

If we wish to derive stellar abundances from dynamic atmo- 
spheres, we have to search for features that are far more sen- 
sitive to changes in the chemical composition than to dynam- 
ical efl'ects. An example would be the ^^CO 3-0 band head. 
As described above, the loop (Fig. [15]) covers a range in the 
colour vs. equivalent width diagram similar to a change in log^ 
by 0.5. This is clearly less than the change from C/O=0.48 to 
C/O=0.25 (Fig.[7]). Therefore we have a good chance of deter- 
mining whether, e.g., the surface abundances of an AGB star 
have already been modified from the first dredge-up composition 
by the third dredge-up (e.g. Lattanzio & Wood 2004). Having 
measurements of diff'erent phases would naturally help us to im- 
prove the accuracy of the result. 

We emphasize that our approach includes some simplifica- 
tions that need to be explored in more detail. We did not test the 
eff'ect of a variation in metallicity. Furthermore, a change in C/0 
may have a more complex eff'ect on dynamic atmospheres than 
known from the hydrostatic case. Finally, a comparison with ob- 
served time series of a mira reveals that some features are in 
common with the model (e.g. the occurrence of loops), while 
other aspects are diff'erent (e.g. the dependency on phase). At 
present, it is not clear if this is due to a general problem in the 
models or a dependence of these characteristics on the param- 
eters of the model (note that RCas is very diff'erent from the 
models we used in our analysis). However, we point out that 
the values found for the equivalent widths in R Cas are in good 
agreement with what we would expect from an extension of the 
detected trends towards redder colours. 

In a forthcoming paper, we will continue our comparison of 
high resolution synthetic spectra from dynamical and hydrostatic 
model atmospheres to identify the most well-suited features for 
measuring stellar parameters and abundances for cool large am- 
plitude variables. We also plan to extend the study to dynamic 
models with diff'erent composition. Finally, we intend to directly 
compare the observed spectra of AGB variables in the globular 
cluster 47 Tuc with the results of hydrostatic and dynamic model 
atmospheres. 
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